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ABSTRACT 



The IEM Company's Continuous Systems Modeling Program 
was used to simulate the longitudinal flight control system 
of the F/A-18 aircraft. The model is intended for use in 
investigations cf aircraft response to flight conditions 
which approach spin or stall and is restricred to the 
automaric flaps up (AFU) flight mode. Program outputs 
includ® stabilator deflection, leading and trailing edge 
flap positions, and cross-axis interconnect signals. 
Various stick forces, motion sensor inputs, and air pressure 
inputs were simulated to produce transient control surface 
responses. These computer generated responses exhibited 
characteristics corresponding to predicted aircraft control 
surface movements. 
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I 



INTRODUCTION 



The Navy has experienced the loss of numerous aircraft 
during recent years due to unintentional departure from 
controlled flight. The increased cost and complexity of 
modern aircraft utilizing fly-by-wire flight control systems 
have placed a renewed emphasis on understanding the 
oerfcrmance character ist ics of these aircraft, especially 
near the limits of the flight envelope. The purpose of this 
thesis is to investigate a method of non-real time computer 
simulation of the longitudinal flight control system of the 
F/A-18 aircraft. Future thesis researchers at the Naval 
Postgraduate School will complete corresponding simulations 
of the lateral and directional control systems. The 
combination of these simulations with an existing 
aerodynamic simulation program will yield a complete 
aircraft stability and control model. The primary purpose of 
the model is to investigate methods of designing control 
augmentation systems which actively inhibit or prevent 
departures from controlled flight. Other uses of this model 
would include the capabilities to test new programmable 
memory configurations, to evaluate new components such as 
ootimal observers, to simulate degraded flight conditions 
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(such as a damaged flight control surf ace ) and to recreate 
flight conditions during post accident investigations. 

The method which was chosen for accomplishment of the 
objectives of this J "h es is was the Continuous Systems 
Modeling Program (CSMP) , developed by the IBM 
Company [Ref* 1]. CSMP is a software package designed to 
simulate dynamic systems described in terms of differential 
equations and block diagrams normally encountered in systems 
theory. CSMP allows programming flexibility through the use 
of thirty-four pre-programmed functional blocks which are 
similar tc FORTRAN subroutines. These blocks provide rapid 
access to mathematical functions, switching functions, 
signal sources, logic functions and FORTRAN functions. Since 
this thesis represented the first attempt at the Naval 
Postgraduate School to accurately model the flight control 
system of a modern, highly augmented tactical aircraft it 
was deemed important to concentrate on the physical systems 
rather than become involved in the complexities of numerical 
a n a 1 y si s . 

Alternatives to CSMP which were considered included 
analog programming and FORTRAN programming. Analog 
programming was not selected because it is less 
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representative of the systems to be modeled and it is less 
accurate than CSMP. The concept of programming directly in 
FORTRAN was carefully considered. Since FORTRAN is the 
source language for CSMP the capabilities of CSMP are a 
subset of the capabilities of FORTRAN itself. Additionally, 
CSMP has restrictions on the number of allowable statements, 
constants, variables and ether parameters. Unlike FORTRAN, 
when functional blocks in CSMP are used, the programmer has 
no direct control of mathematical eperations internal to the 
functions. Tte primary reasons for which CSMP was selected 
were its simplified input statements, output statements and 
urogram control statements which facilitated rapid program 
writing and testing. Additionally, the automatic time and 
amplitude scaling, data formatting, and compatibility with 
graphic display devices which CSMP provides are well suited 
for prototype program development. 
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II . MS THODOLOGY 



A. PROGRAM ORGANIZATION AND NOMENCLATURE 

A detailed description of the flight control laws of the 
F/A-18 aircraft is given in the McDonnell Aircraft Company's 
system design report [Ref. 2]. The computer program 
developed in this thesis is based upon figure 16.1 of this 
report, entitled F/A-18 Longitudinal, Mechanical, CAS, and 
DEL Control Law Mechanization. This figure contains six 
pages of block diagrams depicting generation of longitudinal 
control signals which are valid for all aircraft 
configurations and failure modes. The version of flight 
control program incorporated in the flight control computer 
programmable read only memory (PROM) utilized in this 
simulation is 8.2.1, dated August 31, 1982 . 

A brief discussion of the F/A-18 flight control system 
is necessary to facilitate a discussion of control law 
modeling. All computations of ccntrol laws are accomplished 
independently by four channels of digital computation. 
Primary control in the pitch axis is provided by symmetric 
deflection of the horizontal stabilators. Full span leading 
edge flaps and trailing edge flaps are scheduled to provide 
maximum lift to drag ratio during maneuvering, high angle 
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of attack and cruise configurations. Roll control is 
accomplished by ailerons, differential stabilators ar.d 
differential leading and trailing edge flap deflection. 
Directional control is maintained by dual rudders. The 
thesis computer program which is contained in Appendix A 
simulates the output cf one channel of digital computation 
and calculates the angular oositions of the stabilators, 
leading edge flaps and trailing edge flaps. This simulation 
does not calculate rudder or aileron deflection, however, 
all electrical signals reguired for cross axis interconnects 
are provided. 

The task cf programming the information given in the 
longitudinal control law mechanization schematic was 
simplified by two means. First, the program restrictions and 
assumptions tc conditions cf flight which are discussed in 
part B were applied. As a result sections of figure 16.1 
which apply tc mechanical control laws and spin modes, for 
example, were deleted. This reduced the number of schematic 
blocks to be modeled by approximately one third. The second 
simplification arose through a system of nomenclature in 
which nine control paths were defined in order to limit the 
number of input and output signals for any specific path. 
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When combined, these nine paths fcrm the total longitudinal 
control law mechanization. 

Figure 2.1 is the overall longitudinal control signal 
block diagram which was used in this simulation. It was 
derived by applying all program assumptions or restrictions 
outlined in par* B to figure 16.1 of the system design 
re port. 

Block diagrams depicting the logic development of the 
component paths are included in Appendix A. The nomenclature 
for each control path which is given in Table I is peculiar 
to this simulation program. 

Table II lists nomenclature for control signal groups which 
are common to both the McDonnell schematic and to this 
simulation. Signals with common prefixes are numbered 
consecutively in the feed forward direction. The primary 
feed forward input to the CSMP simulation is pilot stick 
force. Feedback signals include Ditch rate, roll rate, yaw 
rate, angle of attack, normal acceleration and differential 
control surface commands. 

B. PROGRAM RESTRICTIONS AND ASSUMPTIONS 

The task of obtaining a working computer program to meet 
the thesis objectives did not require the simulation of all 
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Figure 2.1. 



Longitudinal Control Laws 
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TABLE I 



Ccntrcl Signal Paths 



PATH 



SIGNAL PREFIX 



Mam path (open Icon) MP 

Pilot smck input path PS 

Pitch rate gyro path PR 

Angle of attack sanscr path AA 

Normal accelerometer NZ 

Forward integrator path FI 

Stabilafor path ST 

leading edae flap path LE 

Trailing edge flap path TE 



TABLE II 



Common Signal Nomenclature 



SIGNAL (OR PREFIX ) 



D ESC RI PTION 



ALPHAS 

ALPHA! 

F 

PK 

PS 

PV 

QC 

RI 



Angle of attack (computed) 

Angle of attack (true) 

Function (air data schedule) 

Air data schedule gain parameter 
Static pressure 
Pitch axis storage location 
Compressible dynamic pressure 
Pressure guotient (=QC/PS) 



possible conditions of flight. Thus, an assumed aircraft 
configuration and flight condition led to many 
simplifications in the model. Each major assumption is 
discussed below. Should future researchers desire to 
construct a mere general mcdel, additional program logic 
paths could be included in this CSMP simulation without 
requiring major revisions to the program. 
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Assum pticns : 

1. The aircraft flight control electronics set (FCES) is 
operating normally. Primary components of this set include 
flight control computers, pitch, roll and yaw rate gyros, 
accelerometer assemblies and flight control panels. The FCES 
contains logic sequences for failure detection and 
corrective control law implementation which were ommitted 
from the thesis computer program, 

2. The aircraft is under inner loop control. In rhis 
mode pilot stick force forms the primary control input. 
Autopilot functions such as heading hold, roll or pitch 
attitude held, and altitude hold are not in operation. 

3. The pitch control augmentation system (CAS) is 
oroviding longitudinal control. BackuD control systems which 
are net in opera-ion include the mechanical backup mode and 
the direct electric link (EEL) mode which provides an open 
loop signal from the pilct stick position sensor to the 
stabilator servoact uator . 

4. The aircraft is operating in the up and away flight 
phase in the auto flaps up (AFU) configuration. This phase 
requires that the flap switch in the cockpit be in the AUTO 
position or that the calibrated airspeed be greater than 243 
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knots. In the AFU mode gain schedules are optimized for 
combat maneuvering characteristics in the low to mid dynamic 
pressure regime. Additionally, a trim integrator controls 
load factor and, at angles of attack above 22 degrees, 
proportional ncse down commands are introduced. 

5. The aain switch is in the normal position. The 
override position of this switch causes flight computers to 
use fixed values for air data schedules and a predetermined 
angle of attack in control law computations. The normal 
nosition of the gain switch allows measured values of air 
data and angle of attack to be used. 

6. The aircraft is operating with weight off the wheels, 
soeed brake in, and with no external stores. 

7. Anti spin functions are not simulated. 
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III. CONTROL LAW MODELING 



A. AIR DATA SCHEDULES 

Air data schedules ars functions of sxatic pressure 
(PS) , compressible dynanic pressure (QC) , and other 
parameters such as normal acceleration, angle of attack or 
condition of external stores. The longitudinal control laws 
contain inputs from 13 different air data schedules which 
are listed in Table III. 



TABLE III 



Air Data Schedules 



FUNCTION DESCRIP TION ( INPUTS ) 

F12 Pitch forward loop integrator aain schedule (RI,PS) 
F22 Fader on supersonic compensation (QC) 

F23 Snail margin on pitch forward loop integrator (AOA) 
F24 Trailing edge flap schedule (AOA, RI) 

F25 Trailing edge flap schedule (QC limit) 

F27 Leading edge flap schedule (AOA,RI) 

E28 Leading edge flap schedule (QC limit) 

F29U Leading edue flap schedule (RI limit) 

E32A Longitudinal forward loop gain schedule (QC) 

F37 NZ limit on AOA feedback (NZA) 

F40 Pitch rate feedback gain schedule (PS,QC,RI) 

F68 Pitch rate feedback gain schedule (QC) 

F107 Longitudinal inertial gain schedule (QC) 



Each schedule controls an output gain for a particular 



purpose. Function 32A, for example, yields a uniform initial 



oitch acceleraxicn in response tc sharp inputs, with gain 



decreasing a- high values of compressible dynamic pressure. 



Other functions such as function 29U are used to determine 
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an upper cr lower limit to an input signal. Simulation of 
these functions in the CSME format required ohat individual 
computer programs be written to test each air data schedule. 
In all cases the final result was given in the flight 
control system design report, usually in graphical form. 
Conversion to CSHP format was done by extraction of data 
ooints from graphs, use of logic flow charts when available 
and direct employment of those mathematical formulae which 
were given. Cnee programmed in CSMP, the tabular and 
graphical output data was compared to data given in the 
design report for the same test parameters. Thus, the output 
of each function was independently verified before the 
function was included in the longitudinal control law. 
si mulation . 

The simulation is written tc take advantage of CSMP’s 
SORT and NCSORT capabilities where appropriate. The NOSCRT 
motion of CSHP is used for conditional logic and branching. 
This statement allows the user flexibility in creating 
sections cf the program in which ordinary FORTRAN rules can 
be used. In NOSORT sections intermediate variables were 
defined, comparisons made and appropriate branching was 
executed. The program statements were returned to SORT 
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12 and 40 could have been combined. This would have 
increased difficulty in program debugging and made areas of 
the program which reguire combinational logic less apparent 
to readers. 

Appendix C contains documentation for the air data 
schedules. Each schedule is included in a complete computer 
program which produces tabular and graphical data to match 
the characteristics of figures in the design report. In 
incorporating these air data schedules into the longitudinal 
control law simulation all schedules except for function 24 
were placed before the main computational body of the 
program. Function 24, the trailing edge flap schedule, 
reguires conditional logic and contains computed angle of 
attack as an input. Since computed angle of attack is 
generated in the normal accelerometer path, function 24 was 
placed immediately before the trailing edge flap path 
computations . 

B. FADERS, A NA ICG TO DIGITAL CONVERTERS, MECHANICAL 

EEEAKODTS 

The purpose of faders is to eliminate large 
discontinuities in gain, permitting gradual change in output 
from old values to new new values at a desired sample rate. 
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Fa de rs are located in Figure 16. 1 at the outputs of air data 
schedules and gain schedules which are dependent upon 
aircraft configuration. Discontinuities may arise as a 
result of a change in physical measurements such as dynamic 
pressure, change in aircraft configuration such as speed 
brake extension, or change in mode of flight. The lower 
limit for signal MP4, for example, changes from a gain of 
-50.0 to 0.0 as a result cf spin entry. Since the thesis 
computer program assumes that electrical signals vary 
smoothly and is restricted to up and away controlled flight, 
faders are modeled in freguency only. 

Analog to digital conversion and frequency matching are 
obtained by a sample and held process. Sampling times are 
generated by the CSMP functional block IMPULS which produces 
a time series cf unit impulse functions with a specified 
start time and period. Since these pulse trains are used in 
several areas cf the program, impulse functions of 20, 40 

and 80 hertz are included immediately following the air data 
schedules. The zero order hold function ZHOLD keeps signal 
gain constant throughout the pulse period. 

A mechanical breakout force cf two pounds is modeled in 
the pilot stick input path. When large stick forces are 
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aoplied, the resultant signal is calculated with a reduction 
in magnitude cf two pounds. This modeling is accomplished in 
CSMP by the deadspace (DEADSP) functional block. Figure 16.1 
of the design report depicts electrical signal "deadbands” 
which are conceptually identical to mechanical breakouts. 
Appendix D cortains a program which displays deadspace 
outputs in tabular and graphical form. 

C. ALIASING FILTERS AND SIGNAL LIMITERS 

The longitudinal control laws cf the F/A-18 include five 
aliasing filters which are modeled as first and second order 
Laplace transforms. The first order transform is of the 
tvpe A/(Bs*1) and is in the pilot stick input path. This lag 
type filter with one real pole is converted to CSMP format 
by the functional block REALPL. A required initial condition 
is the value cf the output signal when time is zero. This 
may be determined arbitrarily by the user, but was set to 
zero for this simulation. 

Second order filters are present in the pilot stick 
input path, pitch rate gyro path, angle of attack sensor 
oath and normal accelerometer path. Each is of the form 
A/(Bs 2 +Cs+D) and represents an underdamoed system. The CSMP 
functional block for complex poles (CMPXPL) is used. 
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Initial conditions are tfce value of the output signal and 
time rate of change of the output signal when time equals 
zero. Natural freguencies varied from 4.34 hertz in the 
pilot stick input path to 33.3 hertz in the angle of attack 
sensor path. The appendix contains a computer program which 
demonstrates tfce CSMP outputs of both first and second order 
aliasing filters. Unit step inputs were introduced to each 
filter to generate transient responses. The first order 
filter produced an exponential rise to steady state with the 
correct time constant. C haracte ristics of the second order 
filters such as rise time, peak time, maximum overshoot, and 
settling time compared favorably with theoretical 
results [Ref. 3]. 

Sianal limiters restrict the maximum or minimum values 
of an output signal. Stafcilacor surface deflection, for 
example, is limited to 10.5 degrees trailing edge down and 
24 degrees trailing edge Up. The CSMP functional block LIMIT 
allows direct speci f ica t ion of lower and upper signal 
li mi ts. 

D. FREQUENCY AV ERAG ERS AND RATE LIMITERS 

A characteristic of tfce F/A-18 flight control system is 
that various signal paths operate at freguencies of 20, 40 
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or 80 hertz. When signals are combined mathematically, the 
inputs are first converted to a common frequency. Normal 
accelerometer path signals are computed at 40 hertz and 
combined with outputs of the forward integrator path, which 
operates at 20 hertz. In this case a 20 to 40 hertz averager 
is present between the integrator path and the summing 
junction connected to the normal accelerometer path. The 
algorithm used for the 20 to 40 hertz averager is based on a 
Drocedure given in the Flight Control Electronic System 
Report [Ref. 4]. The averager was required to generate 
signal values at twice the rate of the incoming pulses by 
linear interpolation between amplitudes of the two previous 
signals at 20 hertz. The formulas used to generate 
intermediate signals were: 

Z 40 = Z40Z 1 + DEL (3.1) 

DEI = ( Z 20 -Z20Z1) / 2.0 (3.2) 

where : 

Z40 = current value of the 40 hertz signal 

Z40Z1 = previous value of the 40 hertz signal 
Z20 = current value of the 20 hertz signal 

Z20Z1 = previous value of the 20 hertz signal. 
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Conditional logic was used to keep the output signal 
equal to the input signal an times when the 20 hertz impuse 
function was equal to one. The algorithm was initialized by 
letting "previous values" be equal to input values when time 
equals zero. 

The presence of numerous frequency averagers in the 
control laws would have reguired the repetition of many 
statement blocks without the use of CSMP program MACROS, 
which are similar to FCETRAN subroutines. Frequency 
averagers in tte thesis computer simulation were included in 
MACROS and placed at the beginning of the program. A MACRO 
may be used several times within a program. Input and output 
variables are given dummy names, yet the MACRO is invoked 
with a unique name which is assigned in a function 
definition statement. A limitation to the use of MACROS is 
that variables which are defined in MACRO structure 
statements are not available for output unless they are 
designated as arguments in the function definition 
statement. Additionally, certain functional blocks such as 
REALPL, CMPXPI and INTGRL cannot be used as arguments or as 
parts of a MACRO structure statement. 



27 



The PROCEDURE function of CS MP was used in each MACRO to 



cause statements to be executed in the order of their 
appearance. Each structure statement of a MACRO will be 
individually sorted unless PROCEDURE is specified. All 
statements contained within a PROCEDURE function are treated 
as a single block which can be moved but not rearranged by 
the CSMP translator. 

FORTRAN subprograms were not used in this simulation 
because the size cf the longitudinal portion of the program 
did not approach the maximum limits of CSMP. Inclusion of 
the lateral and directional control systems in this program 
will require that measures be taken to remain within the 
allowable number of structure statements, NOSORT sections 
and MACROS. Tfce CSMP translator is not used to process 
FORTRAN subprogram statements. Since the number of 
subprogram statements is net counted the overall size of the 
oregram may be increased. The capabilities of the computer 
system library may be utilized by subprograms through the 
use cf the CALL statement. The CALL statement must, however, 
be included in a NOSORT or PROCEDURE section. A method to 
invoke subprograms exists which does not require NOSORT or 
PROCEDURE sections, but it is valid only for two or more 
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output variables and it requires a specific format for 
arguments. A final restriction to the use cf FORTRAN 
subprograms is that certain CSMP functional blocks such as 
ZHOLD, IMPULS, and CMFXPL are not allowable. 

The longitudinal control laws contain rate limiters 
which operate at frequencies of 20 hertz and 80 hertz to 
restrict the speed of leading edge flap movement. The 
algorithm for these limiters compares the value of each 
incoming signal with the value cf the previous output 
signal. The magnitude of the difference between these 
signals is processed by the LIMIT functional block which 
generates the current output signal. The 20 hertz rate 
limiter, for example, allows a maximum change in output 
signal value cf 0.9 degrees during each period of 0.05 
seconds. A listing cf the program which was used to test the 
rate limiter is contained in appendix D. Correct outputs 
were observed for both increasing and decreasing input 
signals. 

E. EIGITAL FILTERS, DIGITAL TO ANALOG CONVERTERS, 

SERVOMECHANISMS 

Three types of digital filters are used in the 
longitudinal control laws. Z filter number P2 is a lead-lag 
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type controller in the pitch rate gyro path which operates 
at 20 hertz. Z filter number P9 , the forward loop 
integrator, operates at 2C hertz and compares the aircraft 
response to the maneuver command. The output signal drives 
the stabilator servo act uatcr to reduce the maneuver error to 
zero. This allows the aircraft to to be automatically kept 
in a hands off condition since the forward integrator 
eliminates unccmiranded normal acceleration. Z filter number 
P8 is a structural notch filter which operates at 80 hertz. 
It attenuates aeroelastic bending which is detected by the 
motion sensors. 

In the simulation each filter was developed in its most 
aeneral form for inclusion in a MACRO. The following 
equation for lead-lag filter number P2 is given in the 
Schematic Design Report [ Bef . 5]. 

PR4 (1 + FK1 1*(1-EK12) )Z - ( PK 1 U 1 ) * ( 1-PK 12) 

FR3 ~ Z - (1-PK12) (3.3) 

It is modeled in the thesis simulation as 



FOUT _ AZ - B 

FIN ' Z - C (3.4) 

where A, B and C are constants. The right shifting and 
linearity properties of the z transform are used to solve 
explicitly for the variable FOUT fRef« 6). 
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FOOT 



A - B (Z- 1) 



FIN 1 .0 - C (Z-1) (3.5) 

Cross multiplication and rearrangement yields: 

FO OT= A* FIN - B*FIN (Z-1) + C*FOUT(Z-1) (3.6) 

which is described in the simulation as: 

FOOI= A*FIN - B *FINZ1 + C*FOUTZ1 (3.7) 

This method, which is termed direct realization 
programming, was also used in the development of the notch 
filter and the forward loop integrator. The equivalent 
Laplace transform for each longitudinal flight control 
filter is listed in the design report. This permitted a 
cress check of z filter performance which is included in 
appendix D. 

A specific method of integration may be specified in the 
terminal portion of a CSME program. In the case of flight 
control simulation the Runge-Kutta Fixed Step Size (RKSFX) 
method was utilized to ensure that integrations would only 
occur at the desired sampling rate. The highest sampling 
frequency in ary axis of the F/A-18 flight control system is 
80 hertz, thus the CSMP integration interval DELT was 
specified as 0.0125 seconds. 
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The variable KEEP is used in CSMP to indicate that the 
end cf a valid integration step has bean reached. KEEP is 
set equal to one when this condition is met. During trial or 
intermediate integration steps KEEP will equal zero. Each 
MACRO contains conditional loqic which allows calculations 
to be performed only when KEEP equals one. 

Conversion cf signals from analog to digital form in the 
F/A-18 occurs as the signal reaches the servomechanism. The 
quantizer functicnal block (QNTZ R) is employed to accomplish 
the analog to digital simulation. The transfer function- 
used by the stabilator and flap servomechanisms are not 
oublished ty the manufacturer of the aircraft. In order for 
the thesis computer program to generate control surface 
positions the response characteristics of a Parker Hannifin 
fly-ty-wire actuation system were incorporated. The selected 
actuators were designed for use with all-digital flight 
controls and are modeled as second order systems. The 
stabilator transfer function is 

X _ 1600 

Ei S 2 + 56S + 1600 (3.8) 

where X is the actuator position in degrees and Si is the 
position command. The transfer function for both leading 
and trailing edge flaps is 
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X 

Ei 



400 



(3.9) 



S a + 28S +4 00 

The advantages of using the second order model instead of a 
first order model are that the faster rise time more closely 
represents the physical actuators and that the natural 
frequency and damping ratio may te indeoen dent ly modified. 

F. PROGRAM TESTING METHODS AND RESULTS 

The thesis computer program was tested on three levels. 
The lowest level involved evaluation of the individual 
signal blocks of figure 16.1 in the system design report. 
Sections A through E of this chapter describe signal block 
modeling techniques and appendicies C through G contain 
testing programs. Tc obtain verification of proper program 
operation it was desired to create input signals of 
realistic value which would produce a time varying output. 
The rate limiters, for example, were tested with an input 
signal which rose exponentially to a limiting value, then 
decreased exponentially . Output signals at the desired 

frequency were observed for incoming signals of positive or 
negative slope and for incoming signals which were within or 
beyond the rate limit. 
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The second level of program testing involved the nine 
signal paths which are listed in Table 1. In most cases a 
step function was used as the input signal. Intermediate 
signals were observed to determine continuity, 
frequency , time constants, and conformance with limiting 

values. 

The highest level of program testing required the 
generation of control surface deflections for specific 
combinaticns of pilot stick force and motion sensor feedback 
signals. Neither aircraft flight data nor the McDonnell 
flight control simulation data were available for direct 
comparison with the outputs of the thesis program. The 
description of control characteristics given in the design 
renort was used to make a qualitative analysis of program 
cper ation . 

One specific condition of flight was selected as a basis 
for comparison of control responses to various combinations 
of input signals. This "fcase condition" of flight was 
defined such that inputs from motion sensors could be 
superimposed on base condition inputs to enable an 
investigation cf the effects of each parameter. The base 
condition for program testing was selected to model an 
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Figure 3.1. Base Condition Stabilator Response 

aircraft operating at 20,000 feet and 250 knots. This fixed 
the values of static and dynamic pressures for each 
simulation and permitted manual verification of the gains 
produced by the air data schedules. Additionally, a step 
function representing six pounds of force in the aft 
(positive) direction on the control stick was applied at 
time 0.0 seconds. All motion sensor inputs were held at zero 
so that their effects could be individually studied. The 
initial deflections of the stabilator, leading edge flaps 
and trailing edge flaps were set tc zero degrees and all 
initial conditions for filters were set to zero. 
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Figure 3.2. Base Condition Leading Edge Flap Response 

Figure 3.1 depicts the move ment of -he stabilator in 
response to the base flight conditions. A transient 
oscillation is produced in the first second which results 
orimarily from the second order filter in the pilot stick 
dynamics path. After this oscillation has decayed the 
stabilator continues to deflect at a nearly constant rate, 
since feedback inputs are supressed. The stabilator reaches 
the limiting value of 24 degrees trailing edge up after 9.6 
c e co n Qo . 
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Figure 3.3. Leading Edge Flap Response to AOA Feedback 

Figure 3.2 is a plot of leading edge flap deflection 
(LEFLA?) versus time for the base flight condition. A steady 
state flap deflection of 4.8 degrees is achieved after 1.5 
seconds. Since LEFLAP is a function only of angle of attack, 
static pressure and dynamic pressure it was desired to 
observe the variation in LEFLAP with angle of attack. A ramp 
type increase in angle of attack sensor input (AA1) was 
superimposed upon the base flight condition beginning at 
time 1.4 seconds. The resulting schedule of leading edge flap 
deflection is shewn in Figure 3.3 and is consistent with the 
functional description given in the design report. 
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Figure 3.4. Base Condition Trailing Bdge Flap Response 



Trailing edge flap deflection was modeled similarly in 
Figure 3.4 for the base flight condition and in Figure 3.5 
for the condition in which AA1 is a ramp function starting 
at time 1.4 seconds. 

To determine the effect of forward pressure on the pilot 
stick a step input of -12. C pounds was superimposed upon the 
base flight condition at time 4.0 seconds, which simulated 
an instant anecus reversal cf stick force. Figure 3.6 shews 
that the direction of stafcilator deflection changed abruptly 
upon introduction of the new stick force, and that the 
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Figure 3.5. Trailing Edge Flap Response to AOA 'Feedback 

magnitude cf the steady state rate of stabilator deflection 
was approximately unchanged. Pilot stick force inputs of 
larger magnitude were simulated, but the results are not 
shown. In these cases the rates cf stabilator deflection 
increased while initial oscillatory behavior exhibited 
characteristics similar tc base condition response. 

The influence of angle of attack feedback on stabilator 
oosition is displayed in Figure 3.7. As described in the 
design report, angles of attack in excess of 22 degrees will 
generate a proportional stabilator command causing the nose 



39 



s 




Figure 3.6. Stabilator Besponse to Stick Force Reversal 

of the aircraft to pi-.ch down. In this simulation the angle 
of attack input function AA1 was set equal to 22 degrees 
plus a ramp type increase of one degree oer second starting 
a - time 4.0 seconds. The resulting stabilator deflection was 
identical to that for the base condition until time 4.0 
seconds due to the action of the -22 degree bias in the AOA 
feedback path. The nose dcwn pitching effect of AOA feedback 
is apparent at all later times. In response to the base 
condition the stabilator had reached its maximum limit of 24 
degrees trailing edge down at time 9.5 seconds. The 
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application of AOA feedback restricted stabilarcr deflection 
to 20.9 degrees at time 9.5 seconds. 



s 




Figure 3.7. Angle of Attack Damping 

Normal acceleration damping is shown in Figure 3.8. A 
ramp increase in normal acceleration equal to 1.0 g's per 
second was superimposed on the base condition at time 4.0 
seconds. The output signal varies smoothly due to the fact 
that normal acceleration path outputs are processed by the 
forward integrator path. The stabilator reverses its 
direction cf movement within one second of the time feedback 
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is introduced. In this case the stabilator reached its limit 
of deflection of 10.5 degrees trailing edge down at time 9.0 
seconds . 




Figure 3.8. Normal Acceleration Damping 

According to the design reocrt, the predominant 
contribution to pitch damping is generated by the pitch rate 
gyro path. This is because the pitch rate signal PR 5 is 
summed with the main path signal MP9 downstream of the 
forward integrator as shown in Figure I. Figure 3.9 depicts 
stabilator position versus time for a ramp increase in pitch 
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rate beginning at time 4.0 seconds. The change in direction 
of stabilator movement is much more rapid than that produced 
by angle of attack or normal acceleration damping. 




Figure 3.9. Pitch Rate Damping 

The thesis computer pregram used a significant portion 
of allowable CSMF program size. Table IV, which is extracted 
from the "translation table" section of CSMP output, lists 
the areas of the program which most closely approached the 
size limits of CSMP. 
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TABLE IV 



Program Size Restrictions 



Parameter 


Current 

Program 


CS MP 
Maxi mum 


MACRO and statement outputs 


206 


600 


Statement input work area 


422 


1900 


Parameters-f unction generators 


43 


400 


History and memory block names 


21 


50 


MACRO statement storace 


85 


125 


SORT sections 


5 


20 


Maximum statements in section 


17 1 


600 



IV. CONCIUSICKS AND RECOMMENDATIONS 

I he computer program developed in this thesis simulates 
the performance cf the longitudinal flight control system of 
the F/A-18 aircraft by generating control surface 
deflections and cross axis electrical signals. The responses 
to AOA, normal acceleration and pitch rate feedback 
correspond to the descriptions given by the aircraft 
manufacturer . 

The use of CSMP simplified the task of system modeling 
bv provision cf pre-programmed functional blocks and a 
flexible format for outputs. The CSMP orogram size 
restrictions do not appear to be a factor which would 
prohibit the addition of the lateral and directional flight 
control systems to this simulation. It is recommended that 
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future expansion cf this program be done using techniques to 
conserve program size. FORTRAN subroutines and computer 
library functions should be utilized due to the limited 
MACRO and NOSORT capabilities of the CSMP translator. 
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AFPENDI X A 



FLIGHT CONTROL SYSTEM COMPONENT PATHS 
PILOT STICK INPUT PATH 
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PITCH RATE GYRO PATH 




47 



ANGLE OF ATTACK SENSOR PATH 




48 



NORMAL ACCELEROMETER PATH 
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, I 

POSITIVE LIMITER 



FI4 



FORWARD INTEGRATOR 
PATH 





ns 



DIGITAL FILTER 
P9 



Fit 

, !! 

VARIABLE LIMITER 







50 



STABILATOR PATH 




51 



LEADING EDGE 
FLAP PATH 




52 



LEADING EDGE FLAP 
PATH (CONT'D) 




53 



TRAILING EDGE FLAP PATH 



AA5 (FROM 
AOA SENSOR 
PATH) 



OICITAL filter 
P12 


1 


112 

f 


20 HERTZ 
SAMPLE /HOLD 


1 


1Z21 

* 


X 






- 1 


1Z21X 

f 


vMiAtte 


LIMITER 
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TRAILING EDGE FLAP PATH 
( CONT *D ) 
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APPENDIX E 



FLIGHT COHTBOL COMPUTES PBOGRAM 



* *********** ***** * *** ********* ********************** ***** 
* * * ** 

* ** LONGITUDINAL FLIGHT CONTROL LAW SIMULATION ** 

* * * 

* ***************** *** ****** ****************************** 
* 

* ******************* **** *** ************************* 

* *****INTEGRATING,LAG,AND notch digital filters** *** 

* * ** * # ******** * *** ****** ** * *************** *** ******* 



MACRO FOUT = ZINT (FIN, KA,KB r IMP,FOUTZ1) 

PROCEDURAL 

IF (I MP . NE . 1.0) GO TO 10 
IF (KEEP. NE. 1.0) GO TO 10 
IF (TIME. EQ. 0.0) GO TO 10 
FOUT=KA*FIN + K3*F0UTZ1 
F0UTZ1 =FCUT 
10 CONTINUE 
ENDMAC 

MACRO FOUT=ZLAG(FIN,KA,KB,KC,IMP ,FINZ1 ,F0UTZ1) 
PROCEDURAL 

IF (IMP .NE. 1.0) GO TO 10 
IF (KEEP. NE. 1 . 0) GC TO 10 
IF(TIME.EQ.O.O) GO TO 10 
FOUT= K A* FIN - KB*FINZ1 + KC*FOUTZ1 
FOUTZ 1 =FO UT 
FIN Z 1 = FI N 
10 CONTINUE 
ENDMAC 

MACRO FOUT=ZNOTCH (FIN ,KA ,KE,KC,KD,KE r IMP,FINZ1 , . . . 
FINZ2 . FOUTZ 1 ,FOUTZ2) 

PROCEDURAL 

IF (IMP.NE. 1.0) GC TO 10 • 

IF(KEFP.NE. 1.0) GO TO 10 
IF (TIME. EC. 0. 0) GC TO 10 

F0UT=KA*FIN+KB*FINZ1+KC*FINZ 2- KD*FOUTZ 1 - KS*FOUTZ2 
FOUTZ2 = FCUTZ1 
FOUTZ 1 = FC UT 
FINZ2 = FINZ1 
FINZ1 = FIN 
10 CONTINUE 
ENDMAC 
* 

* ***************************** 

* *****FRECUENCY AVERAGERS***** 

* ***************************** 



10 



0 ) 

0 ) 



GO 

GO 



MACRO Z40=AV2040 (Z20 ,IMP) 
PROCEDURAL 

IF ( KEEP . NE. 1 
IF (TIME .NE.O . 

Z2 0 2 1 = Z20 
Z40Z 1 = Z 20 
DEL=0. 0 
GO TO 10 

5 IF (IMP. EQ. 1.0) 

Z40=Z40Z1+DEL 
GO TC 15 
Z40=Z20 

DEL= (Z2C-Z20Z1) /2 . 0 
Z20Z 1 = Z 2C 
Z40Z1=Z40 
15 CONTINUE 

20 CONTINUE 



TO 

TO 



20 

5 



GO TC 10 
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END MAC 

MACRO Z80=AV4080(Z40 ,IMP ) 

PROCEDURAL 

IF (KEEP . NE . 1 . 0) GO 10 20 

IF (TIME. NE. 0.0) GO TO 5 

Z40Z 1=Z 40 

Z80Z1=Z40 

DEL=0. 0 

GO TO 10 

5 IF (IMP.EQ. 1. 0) GO TO 10 

Z80 = Z80Z1 + DEL 
GO TO 15 
10 Z80 = Z40 

DEL = ( Z40-Z 40Z 1 ) /2 . 0 
Z40 Z 1 = Z40 
Z80Z 1 = Z8C 
15 CONTINUE 

20 CONTINUE 

ENDMAC 

MACRO Z80=AV20SC (Z20 r IMP) 

PROCEDURAI 

IF (KEEP .NE. 1 . 0) GO TO 20 
IF (TIME. NE. 0.0) GO TO 5 
Z20Z 1 = Z 20 
Z80 Z 1= Z20 
GO TO 10 

5 IF (IMP.EQ. 1.0) GO TC 10 

Z80 = Z80Z 1 + DEL2 
GO TC 15 
10 Z80=Z20 

DEL2= (Z2C-Z20Z 1 )/4 . 0 
15 CONTINUE 

Z20 Z 1 = Z20 
Z80Z 1 = Z80 
20 CONTINUE 

ENDMAC 
* 

* ************** ************ ************************* 

* ***** INITIAL CONDITIONS AND CONSTANTS FOR BASE***** 

* *****flIGHT CONDITION OF 2 0,000 FT AND 250 KTS** *** 

* ********** t** ************************************** 

INITIAL 

CONSTANT PS=972. 49. QC=1 12 .73, RAT E20 = 0. 05, R AT E4 0=0. 025. . . . 

PS 3IC1=0. 0 ,PS3IC 2=0. 0 ,PS4IC1 = 0.0 ,PR2IC1=0 . 0 ,PR2lC2=0.0 ,. . . 
AA2IC2 = 0. 0 -NZ2IC 1=0. 0 ,NZ2IC2 = 0. 0 , NZ 4Z 1 =0 . 0 . N Z5Z 1 = 0 . 0 , . . . 

?K 10=0. 46 47, FI6Z 1=0. 0 ,M Pi 1 Z1 = 0. 0 , MP 1 1 Z2=0 . 0 , MP 1 2Z1 =0 . 0 , . . . 
?8A = 0.690 84,P8E=-0. 9 906 8 , E8C=0. 6 6312, P8D=- 0.99068, ... 

DS 1 = 0.0. DLE 1 = 0.0 , DTE 1=0. 0 ,LE2 Z1 = 0 . 0 , TE2 Z1 =0 . 0 , LEDI C1=0 . 0 ,. . . 
LE DIC2=U . 0 , TEDIC 1=0 . 0 ,TEDIC2=0. 0 ,PR3Z 1=0. 0 .PR4Z 1=0.0.. . . 
RATES 0=0. 0 125,AA2IC 1 = 0. 0, PK9=-1 . 154 3 , MP1 2 Z2=0 . 0 ,P8 E=0 . 3539 

* 

DYNAMIC 

* 

SORT 

* 

# ##****###**♦#*** 

* *****iN]:gTS***** 

# ********** ^** ## * 

PS1=6.0*STEF(0.0) 

PR1 =0. 0 

AA 1 = 22. 0 + RAMP (4 .0) 

NZ1 =0. 0 
HR 1 = 0. 0 
YR 1 =0. 0 

* 

* 

* ************************************ 

* *****FUNCTIONS REQUIRING NOSORT***** 

* ************************************ 
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* 

♦ 

* 

* 



5 

6 

* 

* 

♦ 

* 



********************* 
*****F UNCTION 1 2***** 
********************* 
RI=QC/PS 

F12T1=RI**2*9.625-. 025*RI+ 
F12T2=PS*7. 96 9E-4+0 .84 
F12MAX=LIMIT(1 .0,8. 0.F12T2 
F12T3=LIMIT (1 . 0 ,F 1 2 MAX, FI 2 
F12T5=LIMIT (. 5, 1.35.RI) 
F12T6=F12T5*(0. 00 952 *PS + 4 . 
F12T4=LIMIT (1 . 0 ,8. 0 ,F12T6) 
NCSORT 

IF( .NOT.RI.GT. 0.5) GO TO 5 
F 1 2 = F12T4 
GO TO 6 
CONTINUE 
F 1 2 = F12T3 
CONTINUE 
SORT 

** ****** * * ********** * 

*****j>(j ACTION 4 0***** 
*********** ** ******** 

F4 0T 1 = 3.25-3. 0*LIMIT (0. 75, 
F40T 2=0 . 656 25- . 0013 125*LIM 
F40T3=-0.2617 7+ (9.6 35E-4) * 
PIQ=LIMIT (0.0, 1 800. C, PS) 

F4 0T4 = LI KIT (0 . 0 , 1 0 5 . 0 , QC-3 
F40T 5= (F40T4* (-PIQ* (1.1742 
+ 0.475- (6 .5E-4) *PIQ 
F40T6=F40T3+(F40T 2*11 MIT (0 
F4 0T7=F4CT4*(5. 674 6E-4+ 1. 9 
+ JO. 16923-3. 84615E-5 
F4 0T8= (0. 16 92 3- 3. 86 415E-5* 
+ (F40T4 * ( 1. 67247E-3- 



1.0 

^ 1 ) 

04) + (-PS*0. 00396- 1. 18 ) 



0. 85. RI) 

IT(500. 0.1800.0. PS) 
LIMIT (500. 0, 180 0.0,PS) 

35.0) 

8E-6) + (1.5 238E-3) ) ) . . . 

.75,0.85, RI) ) 

84 1 E- 7 *PIQ ) . . . 

*PIQ) 

PIQ)... . 

9. 29152E-7*PIQ) ) 



NCSORT 

IF ( .NOT .RI.GE. 0. 75) GO TO 20 
F40T9=F40T8*F40T 1 
GC TO 30 



20 


CONTINUE 

F 40T9 = F 40T 8 




30 


CONTINUE 






IF ( .NOT . (QC. GE .440 . C.AND. 
F 4 0T 1 0= F40T 6 
GO TO 50 


RI.GE 


40 


CONTINUE 

F40T 10=F40T7 




50 


CONTINUE 






IF ( .NOT .PIQ. GT .980 . C) GO 
F 40=F40T9 
GO TO 80 


TO 60 


60 


CONTINUE 






IF ( .NOT .PIQ. LE .500 . C) GO 
F40=F40T5 
GC TO 80 


TO 70 


70 


CONTINUE 

F 4 C = F 40 T 10 




30 

* 


CONTINUE 





RI.GE. 0.75) ) GO TO 40 



* ************************** 

* *****SORTED FUNCTIONS***** 

* **** **** ** ********* ******* 

SCRT 

* ************************* 

* *****FUKCTICNS 22,23***** 

* ************************* 

F2 2 = 0. 0167* (-800. 0 + LIMIT (8 00. 0, 900.0 ,QC) ) 
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* 

* 

* 

* 

* 

* 

* 

* 

* 

♦ 



* 

♦ 

* 

* 



* 

* 

* 

* 



❖ 

* 

* 

* 



♦ 

* 

* 

* 



F 23 = 3. 14 35 -0. 14 29* LIMIT (15 .0, 21 .998, ALP HAT) 

$**:}e:*:***4** ****** ****** ** ****************** 

***** F UNCTIONS 25 ,2 7,28 ,29 U, 29L, 32A , 37***** 
**#***.#*#$* ****** ************************** 

F25 = 47. 6 56-0. 05 106* LIMIT (6 00. 0,8 35. 0,QC) 

F27= 1.3 2 8* (ALPHAT + 7 . 858 4- 1 7. 86* LI MIT (0 . 4 4 ,0. 63, R I) ) 

F28=44. 551-0. 04 05 8*LI MI T(260. 0,950. 0 ,QC) 

F29 0 = 87 . 382 5- 76. 25*11 MI T(0 .7, 1. 146, R I) 

F2 9 1=0 . 0 

QKF = LI«IT (200. 0,2 00 0.0, QC) 

F32A=10C. 0/QKF 

F37 = 2.5-0. 5*LIMIT(3.C,5.0, NZA) 

********** *** * * ** < <***** 

*****F UNCTION S 68,1 07***** 

****** * * * * ********* * *** *** 

F6 8 = -0.00 2977* (-4 80 . C+LIMI T ( 250 . 0, 48 0. 0 , QC) ) 

F 10711= ( (-5.71 4E-7) *ABS (QC - 750. 0) ) + ( 8. 4E- 4) 

F 1 0 7 =LI KIT (0. 0, 10000.0, F10 7T1) 

***************** ********* * 

*****IMPULSE FUNCTIONS**** * 
*************************** 

IMP20=IMFULS(0. 0,RATE20) 

IMP4 0=IMPOIS(0. 0, RATE40) 

IMP80=IKFUIS(O. 0, R ATE80) 



*** ******* 
***«*PK F u 
********** 
PK 1 1=1 . 65* 
PK1 2=0. 565 
PK 16A=3. 5* 
PK 1 7 A= F37* 
PK 19A=F 12* 
PK2 1 A=F 1 2* 
PK22A=F 12* 
PK1 6 =ZHC ID 
PK 17=ZHCID 
PK1 9=ZHC ID 
PK2 1 =Z HOLD 
PK22=ZHCID 



************ 

NOTIONS***** 

************ 

F22 

4 

F 32 A 
F32A 

F32A * 3. 5 
F32 A* F68 
F 32 A 

(IMP20, PK 16 A) 
(IMP20, PK 17A) 
(IMP20, PK19A) 
(IMP20,PK21A) 
(IMP20, PK22A) 



* 

* 

* 



******************************** 

*****PHCT STICK INPUT PATH***** 
******************************** 

PS2 = EEACSP (-2. 0,2. 0,ES1) 

PS3 = CMPXPL (PS3IC1 , PS3IC2, 0 .14, 27. 3, PS2 *106. 4 7) 
PS4 = REALPL (PS4IC1 , 7.936 5E- 3,PS3) 

PS5 = ZHC ID (IMP80 ,PS4) 

PS6=PS5 

PS7=PS6* (7.0+ (0 .2* AES (PS6) ) ) 

PS8 = LIM IT (-25. 0 ,50. 0,PS7) 

PV3=FS8 

PV3A=ZHCLD (IMP 20, PV3) 

PS9=FV3*F32A 

****************************** 

*****PITCH RATE GYRC PATH***** 
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* * * #*** * * * * * * * 



♦ ****************************** 

PR2=CMPXFL (PR2IC1, PE2IC2, 0 . 8 9, 78 . 5, P R 1 *6 1 6 2. 25 ) 
PR 3 = ZHO ID (IMP80 ,PR2) 

P2A = 1. O + FK 11* ( 1 .0-PK12) 

P2B= (1.0 + FK11) * (1.0-EK12) 

P2C= 1. O-PK 12 

PR4 = ZLAG (FR3. P2 A, P2E.P2C.I MP80,PR3Z 1 ,PR4Z 1) 
PV4=LIMIT (-80 . 0,120. 0.PR4) 

PR 5 = PV 4 * < ( F68 *F 32 A) + F40) 

PV4 A =Z HOLD (IMP20,PV4) 

** * # ** ###*#** * ********* ****** *** * * 

*****ANGIE CF ATTACK SENSOR PATH***** 

*** ********************************* * 

AA2=CMPXFL (AA2IC1 , AA2IC2, 0 . 74 ,209.0, AA 1 *4 368 1.0) 
AA3=ZHOID (IMP40 ,AA2) 

ALPHAT=0.59*AA3+1 . 9 
ALPH AS = AIPHAT 
AA6=ALPHAT-22. 0 
AA7 = IIMIT (0.0 , 100 0 0.0, A A6) 

PV 1 = PK 1 7*A A7 
PV1 A = ZHCID (IMP2 0, PV 1) 

*********************************** 

*****NORKAI ACCELEROMETER PATH***** 
*********************************** 

NZ2 = CMPXFL (NZ2IC1 , NZ 2IC2, 0 . 8 9 , 200 . 0 , NZ 1 *4 000 0. 0) 
N Z3 = ZHO ID (I MP4 0 ,NZ2) 

N Z A = NZ 3 

NZ4 = NZ3- (RR1**2 )* (6.8529E-6) 

?5A = ( 1.0 + PK9* ( 1 .O-PK 10) ) 

P5B= (1.0 + PK9) * (1.0-EK10) 

P5C= 1 . O-PK 10 

NZ5 = ZLAG (NZ4. P5A, P5 E , P5C, I MP40 , NZ 4Z 1 ,NZ5Z 1) 
PV2=LIMIT (-10. 0 ,10. 0,NZ5) 

NZ7=FV2*FK 16 
PV 2A=Z HOLD (IMP20,PV2) 

********************************* 

*****FORWARD INTEGRATOR PA TH***** 

******** ***** * * ** ****** *** ******* 

FI2=MP 3 

FI3 = II MIT (-10 000. 0.0.0, FI2) 

FI4=LIMIT (0.0, 1 0000 .0,FI2) 

FI5= (FI3*F23) +FI4 

FI6 = ZINT (F 15, 0. 05 , 1 . C,I MP2 0,FI6Z1) 

******************* 

***** MAIN path ***** 

******************* 

MP 1 = PV 3A*P K 22 
MP2= (PV 4 A*FK2 1) -MP 1 
MP3= (P K 1 9* PV2 A ) +PV 1 A+MP2 
MP4 = LIMIT (-50 . 0 ,25. C.FI6) 

MP5 A = A V204 0 (MP4 ,IMP20) 

MP5 = ZHOLD (IMP 40 ,MP5A) 

MP6=NZ7+KP5 
MP7=MP6+PV 1 

MP8 A = M P7 + ( (RR1*YR1) *F107) 

MP8=ZHOLD (IMP40 ,MP8A) 

MP9 = AV4080 (MP8 , IMP4 0) 

MP1 0=MP9+PR5 

MP 1 1=LIMIT (-25. 0,25 .C,MP10) 

MP 1 2=Z NCTCH (MP1 1, P8A,F8B,P8C,P8D,P8E,IMP80,. .. 

MP 11Z 1,MP1 1Z2, MP12Z1, MP12Z2) 

MP 1 3 = MF 1 2- FS9 

************************* 

***** STA eilator path***** 
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* * * * * * ■* oo cn * * * * * 



* **** **** 4* *** * *** ****** 4* 

ST2=IIMIT (-24. 0 , 10. 5,MP13) 

ST3 = ST2-DS 1 

ST 4 = LIMIT (-24. 0 ,10. 5,ST3) 

ST5 = CNTZR (0.0 1 2 5, ST 4) 

ST A D EF=CMPXPL (TEDIC 1 ,TEDIC 2 , 0 . 7, 40 . 0 ,ST5* 1600.0) 



4**4* LEADING 2DG g p LA p p A TH***** 

4444*444444444 4 44 444444444 444444 

LE2=ZINT (ALPHAS .0. 0625, 0. 9 375 ,1 MP80 , LE2Z1) 
LE2T=ZHCLD (IMP20,LE2) 

LE2T1A=IE2T*F27 

LE2T1=ZHCLD(IMP20,LE2T1A) 

LE 2 T 2= LIMIT (0. 0 ,F2 8 , LE2T1 ) 

LE3=LIMIT (C.0,F29U,LE2T2) 

44 4 44 4* 4* 44*4***** ****************** ** 

44*** 18 DEG/SEC RATE LIMIT, 20 HZ***** 

4 * *** * * 4* * 4 4 4 ** *** 44**44*44*444*4***4* 

NO SORT 

IF (TIME. RE. 0. 0) GO TC 85 

LE4Z 1=LE 3 

CONTINUE 

DEL 1 = LE 2-LE4Z 1 

DEL UK* II HIT (-0.9 ,0 .9, DELI ) 

LE4A=LE4Z1+DELLIM 
IF (IMP20. NE. 1.0) GC TO 90 
LE4Z 1=LE4A 
CONTINUE 

LE4 = ZHOID (IMP 20 ,LE 4 A) 

44*4* 4* 4* ********* ******************** * 



SORT 

LE5 = AV2080 (LE4,IMP20) 

DLE2P 1=LE5 + 3. 0 

DLE2 = DEADSP (-DLE2P1 , ELE2P1 , DLE1) 

DLE3=DLE 1 + DLE2 

NOTE: POSITIVE LEADING EDGE FLAP PATH ONLY 

LE 6 = LE5 + DLE3 

LE 7 = II MIT (-3. 0 ,33. 0,LE6) 

4 4 **4 4 4 44 44*444 4 * 4444444444444444444 4 * 

444**18 DEG/SEC RATE LIMIT, 80 HZ***** 

44444 4 *** 4 * 44 ** 44 * 444444 ************** 



NO SORT 

IF (TIME. NE . 0. 0) GO TO 110 
LE 8 Z 1= IE7 
110 CONTINUE 

DEL1=IE7-LE8Z1 

DEL L I M = LIM IT (-0.22 5,0.225, DELI) 

LE 8 A = IE8Z1 +DELLIM 
IF ( IMP80. NE. 1.0) GO TO 120 
LE 8 Z 1= LE 8 A 
120 CONTINUE 

LE 8 = ZHCID (IMP 8 0 ,LE 8 A) 

4 444 44 4 444 4*44*4*4* *4* 4*4*4**4***4 44 4 44* 

SORT 

LE9 = LIMIT (-2. 0 . 33. 0 ,LE 8 ) 

LEI 0=CNTZR (0. 0 1 25, LE9) 

LEFL AP=CMPXPL (LEDIC 1 ,LEDIC 2 , 1 . 4 , 20 . 0 , LE 1 0*40 0. 0) 

* 

* ***************************** 



* 

* 



** * * *NOSCRT FUNCTION 24 ***** 
4 ** 4 * 44 * 44 ***** ** ********* *** 



F24I1 = 22.538-20.51*IIMIT (0 . 27 , 0 . 66 , R I) 

F24L2 = 32.76-3 6 . 0*LI MIT ( 0. 6 6, 0 . 9 1 , RI) 

F24T 1 = LIMIT (0. 0 . 1 0 0 00 .0 , AL PH A T) 

F24T2= ALPH AT- (14. 8 7 69-7.69 23* LI MIT ( 0 . 27 , 0 . 9 1 ,RI) ) 
F24T3 = -2.0*LIMIT(0. 0,10 000 .0,F24T2) 

F24T4=1.4* (F24T 1+F24T3) 

NOSORT 
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* # # # 



10 

15 



IF (.NOT.RI.GT. 0.66) GO TO 10 
F24L=F24I2 
GO TO 15 
CONTINUE 

F24L=F24L1 
CONTINUE 

F24=LIMIT(0.0,F24L,F24T4) 

*********** ******** ****** ******** 

♦****TH AmNG EDGE FLAP PATH***** 

*************** ** ****** ** ******** 

TE2 = ZINT (ALPHAS .0.0 3125,0. 096 875 , IMP 80 ,TE 2Z 1 ) 
TE2T=ZHCLD (IMP20,TE2) 

TE2T1=TE2T*F24 
TE 3 = LIMIT (0.0,F25,TE2T1) 

* NOTE :POSITIVE TRAILING EDGE FLAP PATH ONLY 

TE4=TE3*DTE 1 
TE 5 = II MIT (-8. 0 .45. 0 ,TE4) 

TE6=CNTZR (0.0125,TE5) 

TEFLAE=CMEXPL (T ED IC 1 ,TEDIC 2 , 1 . 4 , 20. 0 , TE6* 40 . 0) 

TERM INAL 
ME TH CD RKSFX 

TIMER FINTIM=10. 0 ,OUTDEL = 0 .01 25, PRD EL= 0. 1 2 5 0 , DE LT= 0. 0 1 25 
PRINT STAEEF, L EF LAP, TEFL AP 
LABEL STADEF IN DEGREES 
LABEL TIME IN SECONDS 
PAGE XYPLOT 

OUTPUT TIME, STADEF 
END 
STOP 
EN DJCE 
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AFPENDI X C 



AIR EATA SC HEDULES 

LABEL FUNCTION 12 
* 

INITIAL 

PARAMETER PS = (2 00 0 . C, 1 00 0 . 0 , 50 0. 0 , 20 0 . 0 ) 

DYNAMIC 

e p RT 

RI = R AMP (0.0) 

F12T1=RI**2*9.625-. 025*RI+ 1.0 

FI 2T2=PS*7. 96 9E-4+0 .84 

FI 2 MAX=IIMIT( 1 .0,8. 0.F12T2) 

F12T3 = LIMIT (1 . 0 ,F 1 2 MAX, FI 2T1) 

F 1 2T 5=LIMIT (. 5, 1. 3 5 ,SI) 

F12T6=F12T5*(0 . 00 95 2*PS+4. 04) + (-?S*0. 00396- 1. 18 ) 

F 1 2T 4=LIMIT (1 . 0 ,8 . 0 ,F 1 2 T6) 

NCSORT 

IF(.NOT.RI.GT.0.5) GO TO 5 
F 1 2 = F12T4 
GO TO 6 

5 CONTINUE 

F 1 2 = F12T3 

6 CONTINUE 

* 

TERM INAL 

TIMER FINTIM =2.0, CUTDEL=0.05, PRDEL=0.05 
PRINT F12T1,F 12T2,F 12T3 ,F1 2T4 ,F1 2T5 , FI 2T6 ,F12MAX ,F12 
PAGE • XYPICT 
OUTPUT RI, Fl 2 
END 
STOP 
ENDJCB 



LABEL FUNCTION 24 
* 

INITIAL 

PARAMETER RI= ( 0 . 27 , 0 .66 , 0. 9 1 ) 

* 

DYNAMIC 

* 

S C RT 

ALP HA = RAMP (0. 0) 

F24I1=22.538-20 .51*LIMIT(0 .27,0.66,81) 

F2 4 L 2 = 3 2. 7 6-36. 0*LI MIT (0. 6 6 , 0 . 9 1 ,RI) 

F2 4T1 = LIMIT (0 . 0, 10 000.0, ALPHA) 

F24T2=ALPHA-( 14 .87 69-7. 692 3*LIMIT(0. 27, 0. 91, RI) ) 
F24T3=-2. 0*LIMI T (0. 0,10000 .0, F24T2) 

F24 T4= 1 .4* (F24T 1+F2 4T3) 

* 

NOSORT 

IF (. NOT. RI.GT. 0.66) GO TO 10 
F24L=F24L2 
GO TO 15 
10 CONTINUE 

F24L = F24L1 
15 CONTINUE 

F2 4 = II MIT (0.0 ,F24L, F24T4) 

* 

TERM INAL 

TIMER FINTIM=30.0,OUTDEL=0 . 5 , PRDEL=0 . 5 

PRINT F24L 1 ,F 2 4L2 , F 24T1 ,F2 4T2 ,F2 4T3 , F24T4 ,F2 4L,F 24 



63 



PAGE X YE ICT 
OUTPUT ALPHA, F24 
END 
STOP 
ENDJCE 



LABEL FUNCTION 29 
* 

INITIAL 

PARAMETER PS= (1 250. 0) 

* 

DYNAMIC 

RI = R AMP <0. 0) 

* 

F29U = 87. 38 25-7 6 .25* LIMIT (0 . 7 , 1 . 1 46 , R I) 
F29L=0 . 0 

* 

TERMINAL 

TIMER FINT IM= 1 . 5, PRDEL=0.05, OUTDEL=0.05 
PRINT RI,F29U,F29L 
PAGE XYPLOT 
OUTPUT EI,F29U,F29L 
END 
STOP 
EN DJCE 



LABEL FUNCTION 40 
* 

INITIAL 

PARAMETER QC= ( 1 70 0 . C , 44 0. 0 , 33 5 . 0) 

* 

DYNAMIC 
S C ST 

?S= 0. 0001 + RAMP (O.C) 

RI=QC/FS 

F40T 1 = 3. 25-3. 0* LI MIT (0. 75, 0.85. RI) 

F4 0T2=0.65 625-. 00 1 3 125*LIM IT ( 500 . 0, 1 800 . 0 , PS ) 
F40T3=-0.26177+(9.635E-4) * LIMIT ( 500 . 0, 1 80 6 . 0 ,PS> 
PIQ = LIMIT (0.0, 1 800. C. PS) 

F4 0 T 4 = LI MI T (0 . 0 , 1 0 5 . 0, QC-3 35 . 0) 

F40T5= (F40T4* (-PIQ* (1.17428E-6) + ( 1 . 5 238 E-3) ) ) . . . 
+ 0.475- (6 .5E-4) *PIQ 

F40T6=F40T3+(F40T2*LIMIT (0 . 75 ,0 . 85 , RI) ) 

F40T7 = F40T4*(5. 67 46E-4+ 1 . 9 84 1 E-7*PIQ) . . . 

+ (0. 16923-3. 8 4615 E-5*PIQ) 

F4 0T8= (0. 16 923- 3. 86415E-5* PIO) . . . 

+ (F4 0T4*( 1. 6 7 247E-3-9.29152E-7*PIQ) ) 

* 

NCSORT 

IF ( .NOT. RI.GE . 0.75) GO TO 20 
F40T9 = F40T8*F40T 1 
GO TO 30 

20 CONTINUE 

F 40T S=F 40T 8 
30 CONTINUE 

IF ( .NOT . (QC. GE .440 . C. AND. RI .GE. 0.75) ) GO TO 40 
F40T10=F40T6 
GO TO 50 

40 CONTINUE 

F40T 10=F40T 7 
50 CONTINUE 

IF (. NOT. PIO.GT. 980.0) GO TO 60 
F40=F40T9 
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60 



GO TO 80 
CONTINUE 

IF (.NOT.PIO.LE.500 .C) GO TO 70 
F4G=F40T5 
GO TO 80 
70 CONTINUE 

F 4 0 = F40T 10 
80 CONTINUE 

* 

TERMINAL 

TIMER FI NT I M= 2 000. 0 ,PRDEL= 20. 0, OUTDEL = 20 . 0 
PRINT PS,QC,RI , F4 0T7,F4 0T9 ,F40T1 0,F4 0 
PAGE XYPIOT 
OUTPUT PS,F40 
END 
STOP 
END JCE 



LA 

* 

* 

IN 
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* 



* 

TE 



EN 



BEL GENERATED FUNCTION 

THIS PROGRAM GENERATES QC FUNCTIONS 22, 25, 28, 32A, 
ITIAL 

PARAMETER FS=1000.0 
NAMIC 

QC=R AMP (0. 0) 

F2 2=0. 0167* (-80 0. 0 + LIMIT (8 00. 0, 900. 0 ,QC) ) 

F 25 = 4 7. 6 56-0. 05 106* LIMIT (6 00. 0 ,835. 0 ,QC) 

F2 8 = 4 4. 551-0. 0 4 05 8* LIMIT (2 60. 0,9 50. 0 ,QC) 

QKF = LIMIT (200. 0 ,2 0 0 C*. 0 , QC) 

F32A=100.0 / QKF 

F6 8 = -0. 002977* (-48 0 .0+LIMIT (260. 0,4 8 0. 0,QC) ) 

RM INAL 

TIMER FINTIM=2500. 0, OUTDEL=50.0. PRDEL=50.0 
PRINT QC ,F22, F25, F28,QKF,F32A,F68 
PAGE XYPLOT 
OUTPUT QC , F22 
OUTPUT QC,F25 
OUTPUT QC , F28 
OUTPUT QC , F 32A 
OUTPUT QC , F68 
END 
STOP 
DJOE 



LABEL GENERATED FUNCTION 
* 

IN ITIAL 

PARAMETER RI= (0 . 4 4 , 0 . 63 ) 

* 

DYNAMIC 

ALPH A= RAMP (0.0) 

NZ A=ALPHA 
ALPH AT= AIPHA 

* 

F23 = 3. 14 35-0. 14 29* LIMIT (15 .0, 21 .998, ALPH AT) 

* 



AND 68 
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F27=1. 328* (ALPHA+7. 858U-17 .86*LIMIT(0.44,0.63,RI) ) 
F3 7 = 2. 5-0. 5*L I MIT (3. 0,5.0, NZ A) 

ERM INAL 

TIMER FINTIM=30 .0 , PRDEL= 1.0, OUTDEL=1.0 
PRINT F23,F27,F37 
PAGE XYEIOT 
OUTPUT ALP HAT , F 23 
OUTPUT ALPHA, F27 
OUTPUT NZA ,F37 
END 
STOP 
SNDJCE 



66 



APPENDI X D 

COMPUTER PROGRAMS FOR SIGNAL BLOCK TESTING 



EEADSPACE FUNCTION AND ALIASING FILTERS 



INITIAL 

CONSTANT FS2C1=2.0 

* 

DYNAMIC 

PS1= -5.0 + RAMP (0.0) 

PS2=D E ADSP (-PS2C1,PS2C1,PS1) 

PS3=C M PXPI (0.0, 0.0,0. 14,2 7. 3 , PS 2* 106. 4 7) 
PS4 = 3E ALPL (0 .0, 7. 9365E-3,PS3 ) 

* 



TERMINAL 
TIMER 
LABEL 
PRINT 
OUTPU 
PAGE 
END 
STOP 
EN DJCE 



FI NTIM= 10 . 0 , OUTCEL=0 . 25, PRDEL=0.25 
STICK DEADSPACE PLOT 
P S 1 , PS 2 

T TIME, PS 1 ,PS2, FS3, PS4 
XYPLCT 
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FREQUENCY AVERAGERS 



MACRO Z80 = AV2080 (Z20 , IMP) 
PROCEDURAL 

IF (TIME .NE.O . 0) GO TO 52 
Z20 2 1 = Z20 
Z80Z 1 = Z 20 
GO TO 54 

52 IF (IMP. EQ. 1.0) GO TO 54 

Z80=Z80Z1+DEL2 
GO TO 56 
54 Z80=Z20 

DEL 2= (Z2C-Z20Z 1) /4 . C 
56 CONTINUE 

Z 20 Z 1=Z 20 
Z80 Z 1 = Z80 



i 



ENDM AC 
INITIAL 

CONSTANT RATE=0.05 
DYNAMIC 

Y 1 = RAMP (0.0) 

Y1 A = EXP ( Y 1 ) 

Y 2= IHPULS (0.0, RATE) 

Y3= ZHC1D (Y2, Y1 A) 

Y4= AV2080 ( Y 3 , Y 2) 

TERMINAL 

METHOD RKSFX 

TIMER FI NT I M= 1 . 0. OUTDEL=0. 0 1 2 5 , PRDEL =0 . 0 1 25 , DELT = 0 . 0 1 2 5 
LABEL 20 TO 80 AVERAGER 
PRINT Y2,Y3,Y4 
END 
STOP 
EN DJCE 

MACRO Z40 = AV2040 (Z20 ,IMP) 

PROCEDURAL 

IF (TIME. NE.O. 0) GO TO 5 

Z20Z 1 = Z 20 

Z40 Z 1 = Z20 

DEL=0. 0 

GO TO 10 

5 IF (IMP.EQ. 1. 0) GO TO 10 

Z4 0= Z4 0 Z 1 + DEL 
GO TO 15 
10 Z40=Z20 

DEL = (Z20-Z 20Z 1) /2. 0 
Z20Z1=Z20 
Z40Z 1 = Z40 
15 CONTINUE 

ENDM AC 
INITIAL 

CONSTANT RATE=0.05 
DYNAMIC 

Y 1 = RAMP (0.0) 

Y1 A = EXP ( Y 1 ) 

Y 2= IMPULS (0. 0 , RATE) 

Y 3= ZHOIE (Y2 r Y 1 A) 

Y4 = A V2040 (Y3, Y2) 

TERMINAL 

METHOD RKSFX 

TIMER FI NT I M= 1 . 0,OUTDEL=0. 025 r PR DEL = 0 . 0 25 r DELT=0 .0125 
LABEL 20 TC 40 AVERAGER 
PRINT Y1,Y1A,Y2,Y3,Y« 

PAGE XYPLOT 
OUTPUT TIME,Y1 A ,Y3 , Y4 
END 
STOP 
ENDJOE 
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BATE LIMITERS 



IN ITIAL 
DYNAMIC 

Y 1 =R AMP (0.0) 

Y2 = E XP ( Y 1 ) 

Y3=L IM IT (C.0, 3.0, Y2) 

Y4 = (-1 .0 + EXP (-TIME* 2.0) ) *STEP (2.0) 

Y5 = Y 3 + Y4 

IMP2 0=1 MEOIS (0 . 0,0 . 05) 

IMP4 0=IMFULS( 0. 0, 0. 025) 

IMP80=IMPULS(0.0,0. 0125) 

Y20 = ZHOLD (IMP20 ,Y5) 

Y40= 2HCID (IMP40 ,Y5) 

Y80 = ZHOID (IMP80 ,Y5) 

NOSORT 

IF(TIME.NE.O. 0) GO TO 80 
LES Z 1 = Y20 
80 CONTINUE 

DEL1=Y20-LE8Z1 
DELLIK=LIMIT(- 0.1 ,0.1, DELI ) 

LE8=IE8Z1+DELLIM 
IF ( IMP 20. NE. 1 . 0) GO TO 90 
LE8Z1=LE8 
90 CONTINUE 

TERMINAL 

METHOD RKSFX 

TIMER FINTIM=5. 0,OUTDEL=0. 01 2 5, PRDEL =0 . 0 1 25, DSLT =0. 0 1 2 5 
PRINT Y3,Y4,Y5, DELI , DEL LIM ,LE8Z1 ,LE8 ,Y20 , IMP 20 
PAGE XYEICT 
OUTPUT TIM E ,Y 20 ,LE 8 
OUTPUT TIME ,Y20 ,LE8Z1 
END 
STOP 
EN DJCB 
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DIGITAL FILTERS 



MACRO FOUT=ZINT (FIN , K A, KB ,IMP ,FO UTZ 1) 

FBOCEDURAI 

IF (IMP .NE. 1.0) GO TO 10 
IF (KEEP. NE. 1 .0) GC TO 10 
IF (TIME. EO. 0.0) GO TO 10 
FCUT = K A* Fl S + KB*FOUTZ1 
FOUTZ 1 = FO UT 
10 CONTINUE 
ENDMAC 

MACRO FOUT=ZLAG (FIN, KA, K B , KC, IMP , FINZ 1 , FOUTZ 1) 

PROCEDURAL 

IF (IMP .NE. 1.0) GC TO 10 
IF(KEEP.NE. 1.0) GO TO 10 
IF (TIME. EQ.0.0) GC TO 10 
FOUT=KA*FIN - KB*FINZ1 + KC*FOUTZ1 
FOUTZ1 =FOUT 
FINZ 1 = FIN 
10 CONTINUE 
EN DM AC 

MACRO FOUT = Z NOTCH (FI N .K A , K E, KC, K C , K E , IMP , . . . 

FINZ1 ,FINZ2,FOUTZ1 ,FO UTZ2) 

FBOCEDURAI 

IF (IMP .NE. 1.0) GO TO 10 
IF (KEEP. NE. 1 .0) GC TO 10 
3F(TIME. EQ.0.0) GO TO 10 

FOUT=KA*FIN+KB*FINZ 1 +KC*FINZ 2-K D* FOUTZ 1-KE*FOUTZ2 
FOUTZ 2 = FOUTZ1 
FOUTZ1 = FCUT 
FINZ 2 = FINZ 1 
FINZ 1 = FIN 
10 CONTINUE 
ENDMAC 
IN ITIAL 

PA RAM INTN=.0125 .INTD = 1 . 0 , AA1 1 A 1 =0 . 0 , A A1 1 A =0 . 0 
PARAM PK9=-1.15uS, PK 10 = 0 . 464 7, A A 1 OB 1 = 0 . 0 , AA 1 1 B 1=0 .0 
PAR AM NA=. 6908 4. NB=-. 99 06 8 , NC=. 6 6 3 1 2 , N D=- . 99 068, NE = . 35 396 
PARAM AA IOC 1 = 0.0, AA 1 0C2 = 0 . C , A A 1 1 C 1 =0 . 0 , AA 1 1C2 = 0.0 
LA = (1 . + PK9* (1 .-FK10) ) 

LB= { 1.+PK9)* ( 1.-PK10) 

LC = 1 . -PK 1 0 
DYNAMIC 

A A 9 A = STEF (0.0) 

IM P=IMPULS (0. 0,0.01 25) 

AA 1 0 A = ZHCLD (IMF ,AA9 A) 

AA11A = ZINT (AA10A,INTN,INTD,IMP , AA 1 1 A 1) 

AA 1 1 E=ZL AG (AA10A , LA , LB, LC , IMP , AA 10B1, AA1 IB 1) 

AA11C = Z NOTCH (AA10A,NA,NB,NC,ND,NE,IMP,AA10C1,... 
AA10C2,AA11C1,AA1 1C2) 

TERMINAL 
METHOD RKSFX 

TIMER FINTIM= 1. 5 ,OUTDEL = 0. 012 5,D ELT = 0. 01 25 

PRINT AA10A,AA11 A,AA1 IB, AA11C 

OUTPUT AA 10A (0 .0, 1. 5) ,AA 1 1 A (0 . 0, 1 . 5) , . . . 

AA11B (0.0,1. 5) ,AA1 10(0.0, 1.5) 

END 
STOP 
EN DJOE 
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